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High energy x-ray diffraction is used to investigate the bulk oxygen ordering properties of 
YBa2Cu306+a:. Superstructures of Cu-O chains aligned along the b axis and ordered with peri- 
odicity ma, along the a axis have been observed. For x < 0.62 the only observed superstructure is 
ortho-II with m — 2. At room temperature we find ortho-Ill (m = 3) for 0.72 < x < 0.82, ortho-V 
(m — 5) in a mixed state with ortho-II at x = 0.62, and ortho-VIII (m = 8) at X — 0.67. Ortho-II is 
a 3D ordered structural phase, the remaining ones are essentially 2D. None of the superstructures 
develops long range ordering. Studies of the ortho-II ordering properties in samples prepared with 
x — 0.5 but by different methods show that finite size domains with internal thermodynamic equi- 
librium are formed. The crystal perfection as well as the thermal annealing history determine the 
domain size. The temperature dependence of the observed superstructure ordering is investigated 
explicitly and a structural phase diagram is presented. 



I. INTRODUCTION 

It is now well-established that YBa2Cu306+a; (YBCO) 
is antiferromagnetic (AF) in the tetragonal structural 
phase for < x < 0.35, and becomes superconducting 
at low temperatures in the weakly distorted orthorhom- 
bic phase for 0.35 < x < 1. By high temperature ther- 
mal treatment in a suitable oxygen pressure it is possijJj 
to vary the oxygen composition in a continuous way 
that changes the electronic properties from an AF insu- 
lator via an underdoped to an optimally doped high-T c 
superconductor with T c = 93 K, and a slightly overdoped 
material for x > 0.93. YBCO has therefore become a ma- 
jor model system for basic studies of high-T c supercon- 
ductivity and it is a leading candidate for technological 
applications. 

Structural refinement of neutron diffraction data have 
shown that the variable amount of oxygen resides in the 
CuOz basal planeB. A strong tendency towards for- 
mation and alignment of Cu-0 chains gives rise to the 
orthorhombic distortion below the temperature depen- 
dent transition line between the tetragonal and the or- 
thorhombic phase. The importance of oxygen order- 
ing for the superconducting properties has been veri- 
fied directly from experimental studies where crystals 
are quenched from the tetragonal disordered into the or- 
thorhombic ordered phase. Here it is found that T c of 
quenched YBCO is reduced compared to the equilibrium 
value and increases with tinae|-,when the sample is an- 
nealed at room temperatureaij. Equally, it has been 
observed that the oxygen ordering of quenched YBCO 
crystals increases with timea. The consensus is therefore 
that the Cu-0 chain ordering in the CuOa; basal plane 
controls the charge transfer leading to superconductivity 
in the Cu02 planes. However, in spite of the very large 
number of experimental and theoretical model studies 
there is still no finite microscopic understanding of how 
the charge transfer is controlled by the Cu-0 chain length 
and superstructure ordering. Thus, it is not settled how 
the electronic states formed in the Cu-0 chains hybridize 
with the electronic structure in the Cu02 planes, give rise 
to the charge transfer and contribute to the anisotropy 
of the superconducting properties. A likely explanation 
is that the available structural information is not suffi- 
ciently detailed and unambiguous because the oxygen dif- 
fusion kinetics is too slow at the temperatures where the 
superstructures become stable. It has therefore not been 
established which superstructures are actually formed as 
function of oxygen composition x, impurity level and 
thermal treatment, and how they influence the electronic 
states. 

The orthorhombic phase found below the tetragonal 
to orthorhombic phase transition has the basic ortho-I 
structure, and it is formed by Cu-0 chains that align 
along the b axis with oxygen on the so-called 0(1) site, 
whereas the sites on the a axis (called 0(5)) are essen- 
tially empty. Ortho-I is a 3D long range ordered struc- 



ture, but in commonly prepared crystals true long range 
order is prevented by the formation of twin-domains with 
domain size ranging from a few hundred Angstr0m to 
macroscopic size. Clearly, there is disorder in the ortho- 
I chain structure for compositions x < 1.0. Therefore, 
in thermodynamic equilibrium ordered superstructures 
must be formed for T — > 0. At x = 0.5 an ideal 3D 
ortho-II superstructure may in principle be formed in- 
side the ortho-I twin domains with perfect Cu-0 chains 
on every second b axis and the remaining ones are empty, 
i.e. they contain only Cu. 3D ordering with the Cu-0 
chains stacked on top of one another along the c axis has 
been observed but only with finite size ordering in all 
three crystallographic directions. 

Electron microscopy techniques have had a leading role 
in establishing the superstructures of YBCOoila. How- 
ever, the need for confirmation by bulk structural tech- 
niques is generally recognized, because electron beam 
heating of thin crystals may change the mobile oxygen 
content x and generate transient non-equilibrium sur- 
face structures. Also, it is difficult to obtain quantita- 
tive details about the finite size ordering properties and 
their temperature dependence by these techniques. The 
observed superstructures include the Cu-0 chain type 
of ordering as ortho-II as well as more complex order- 
ing sequences of essentially full (Cu-O) and empty (Cu) 
chains with periodicity ma along the a axis and corre- 
sponding superstructure reflections at modulation scat- 
tering vectors: Q = (nh m 0), where h m = 1/m, and 
n < m are integers, and the coordinates refer to the re- 
ciprocal lattice vectors. Superstructure reflections with 
rn = 2,3,4,5 and 8, which we shall call ortho-II, ortho- 
Ill, ortho-IV, ortho-V and ortho-VIII, respectively, have 
been observed experimentally. Ideally, these superstruc- 
tures may be symbolized by their sequence of full (1) Cu- 
O and empty (0) Cu chains. Thus, ortho-II is simply (10) 
with x = 1/2, and the ortho-Ill sequence is (110) with 
composition x — 2/3. Ortho-IV has an ordered sequence 
of (1110) and composition x = 3/4, and ortho-V is a se- 
quential ordering of ortho-II and ortho-Ill, i.e. (10110), 
with composition x = 3/5. Ortho-VIII is ortho-V com- 
bined with ortho-Ill to the sequence (10110110) and ideal 
composition x — 5/8. In principle similar structures with 
full and empty chains interchanged may be stable, but 
they have not been observed experimentally. ___ 

Superstructures with, unit cells 2y/2a x 2\/2a x di3li3'Eil 
and y/2a x 2y/2a x di£l, the so-called herringbone type, 
have been reported. However, as we shall discuss in Sec- 
tion f\1 these superstructures are most likely not from 
oxygen ordering. The ortho-II and ortho-Ill superstruc- 
tures -have-, been verified.-as— itvulk structural phases by 
x-rayll[L3~E£l and neutrorOErEil diffraction techniques. 
The first obser.jzation by x-ray diffraction was made b« 
Fleming et alx3 for ortho-II and by Plakhty et alEl 
for ortho-Ill. Similarly, the first .neutron diffraction data 
were presented, by Zeiske et q/BH f p r ortho-II, and by 
Plakhty et alx3 and Schleger et al£3 for ortho-Ill. Anal- 



ysis of structure factors obtained from a combination 
of neutron and x-ray diffraction data has unequivocally 
shown that the ortho-II and ortho-Ill supersti 
suit from oxygen ordering in Cu-0 chains^ 
However, relaxation of cations associated with the oxy- 
gen chain ordering contributes significantly to the super- 
structure intensities. In particular the barium displace- 
ment has a strong influence on the x-ray diffraction in- 
tensity. The emplacements .show only minor variation 
with oxygen compositionE3E3. Also, in the ortho-II phase 
there was found no significant change in the displace- 
ments as function of temperature!! 8 !. Compiling previous 
room temperature data from x-ray and neutron diffrac- 
tion studies, we find that the ortho-II superstructure has 
been observed for oxygen compositions 0.35 < x < 0.7 
and the ortho-Ill superstructure for 0.7 < x < 0.77. As 
we shall present below, we have observed clear indica- 
tions of bulk phase ortho-V correlations for x = 0.62 
and ortho-VIII for x = 0.67, but we have found no ev- 
idence for the ortho-IV supestructure. Only a few of 
the previous, studies have been carried out above room 
temperatureQ'E3'Ejn3. 

From previous neutron and hard x-ray diffraction stud- 
ies it has been inferred that the finite size ordering of 
ortho-II results from formation of anti-phase domains in- 
side the ortho-I twin domains-below an ordering temper- 
ature of Ton = 125(5)°CErE3. Anisotropic superstruc- 
ture reflections with a Lorcntzian-squared line shape were 
found in the ordered state! as expected from Porod's law: 



S(q) oc 



for finite size ordered domains with sharp 



boundaries in spatial dimension d. Studies of the order- 
ing kinetics following a quench in temperature from the 
ortho-I into the ortho-II phase have shown a time depen- 
dent domain growth that is algebraic at early times and 
logarithmic at late times. The characteristic time of the 
growth process is activated with an activation energy of 
1.4 eV. At 70 °C it is some days and it extrapolates to 
several years at room temperature. On this basis it was 
suggested that the finite size ordering may result from 
pinning of the domain walls by impurities or defects, but 
recent computer simulations have shown that intrinsic 
slowing down due to the large effective activation en- 
ergy fo*-, movement of long Cu-0 chains may contribute 
as wellEa. 

In the present paper we report on experimental stud- 
ies of the oxygen ordering in YBCO covering the oxygen 
compositions 0.35 < x < 0.87 and temperatures up to 
250 °C, by diffraction of high energy synchrotron radia- 
tion (~ 100 keV). Thus, our studies do not include the 
region at and above the optimal doping level x > 0.93 
where Kaldis et al. have observed structural anomalies 
(see e.g. Ref. p4T). The high energy x-ray diffraction tech- 
nique combines the high penetration power of neutrons 
with the high momentum space resolution. The penetra- 
tion depth of 100 keV x-rays in YBCO is of the order of 
1 mm. This assures that we probe the bulk properties of 
the samples and are insensitive to oxygen diffusion in and 



out of the surface, and studies in sample environments 
with varying temperatures and controlled atmospheres 
are easily accessible. Finally, the synchrotron intensity 
is so high that scattering signals down to a factor of 10 8 
smaller than the fundamental Bragg peaks can be re- 
solved, and the kinetics of the ordering can be studied 
with a time resolution of 1 second. From our studies we 
present temperature scans of the structure factors of the 
superstructures, determine their phase boundaries and 
the nature of the ordering. Firstly, we show that the 
superstructures including ortho-V and ortho-VIII, that 
Beyers et a/.t3 have observed by electron microscopy, 
represent bulk structural phases. Secondly, we present 
extensive studies of the ortho-II superstructure ordering 
in crystals of different quality and thermal treatment. 
Finally, we use the present structural data jointly with 
data compiled from previous studies to establish a struc- 
tural phase diagram that includes the oxygen superstruc- 
tures and the tetragonal to orthorhombic (ortho-I) tran- 
sition temperature, Tt-oi, obtained by neutron powder 
diffraction. We also review and compare with structural 
findings by other groups and discuss our results in rela- 
tion to structural model studies. 

The layout of the paper is as follows: In Section || 
we supply information about the crystal growth and ox- 
idation of the sam ple ( II A ), details about the experi- 
mental setup (II B), and the data analysis (III Cf ). The 
experimental results are presented in Section |III| . First 
we account for th e results for the ortho-II superstructure 
formation ( III A ). This includes the dependence of the 
ortho-II correlation length on crystal quality and ther- 
mal treatment of the sample, the phase transition into 
the ortho-I phase, and the stability range of the compo- 
sition x. Then we describe the ordering properties and 
the stability ra nge of the ortho -Ill superstructure ( IIIBj ). 
In Subsections 1IIC and HID we present the properties 
of the ortho-V and ortho-VIII superstructure ordering, 
respectively, and a structural phase di agram of the oxy- 
gen ordering is presented in Section HIE . In Section 
IV we discuss our experimental structural results in rela- 
tion to other structural findings and their importance for 
charge transfer, and to theoretical model descriptions. A 
concluding summary is given in Section |V|. 



II. EXPERIMENTAL DETAILS 

A. Sample Preparation 

The single crystals used to study the different super- 
structure phases and establish the phase diagram were 
grown in YSZ (ytteia-stabilized zirconia) crucibles by a 
flux growth methodt3 using chemicals of 99.999 % purity 
for Y 2 3 and CuO, and 99.997 % for BaC0 3 . The im- 
purity level of the crystals has been analyzed by ICP-MS 
(Inductively Coupled Plasma Mass Spectroscopy). The 
Zr content of the crystals was found to be less than 10 



ppm. by weight. The major impurities were Al, Fe and 
Zn, the sum of which amounts to less than 0.2 % atom 
per unit cell. When optimally doped (x — 0.93) these 
crystals have T c = 93.2 K and the width of the 10 % - 90 
% diamagnetic response is AT C — 0.3 K. 

Of the four crystals used to study the ortho-II order- 
ing properties as function of sample purity, crystals #1, 
#3 and #4 were grown by the flux growth technique de- 
scribed in Ref. B6L and crystal #2 as described above. 
Crystal #1 was used to study the ortho-II ordering prop- 
erties when exposed to six different annealing methods. 

YSZ crucibles, and chemicals of purity better than 
99.99 % and 99.9 % were used for crystals #1 and #3, re- 
spectively. Crystal #4 was grown in a SnC>2 crucible with 
chemicals of purity better than 99.99 %. The impurity 
level of crystals #1, #3 and #4 has not been determined 
directly, but the flux from the crucibles used to grow crys- 
tals #1 and #4 has been analyzed. Excessive amounts 
of Zr and Hf (up to 15000 ppm for Zr and 550 ppm for 
Hf) were found, but these elements are known to have 
a very low solubility in YBCO. Major impurity compo- 
nents were Al with concentrations 200 ppm and 335 ppm 
in the flux from crystal #1 and crystal #4, respectively, 
and similarly for Eu: 115 ppm and 110 ppm. The super- 
conducting transition temperatures and the widths of the 
transitions have been determined in the fully oxygenated 
state (x = 0.99). They are: T c = 91.0 K, AT C < 1.0K 
for crystal #1, T c = 92.0 K, AT C = 1.5 K for crystal #3, 
and T c = 91.5 K, AT C < 0.5 K for crystal #4. The lower 
T c and larger AT C values of these crystals are not neces- 
sarily a consequence of bad crystal quality but rather a 
result of overdoping. Thus crystal #1 has T c = 93.5 K 
and AT C < 0.5 K when optimally doped. For all crystals 
prepared by high purity chemicals (better than 99.99 %) 
it is likely that the impurities come from minority compo- 
nents of the crucible material or from the furnace walls. 
All the crystals are plate-like with thicknesses of \ to 1 
mm, flat dimensions of 1.5 to 3 mm and weights ranging 
from 10 to 70 mg. 

" 31 s 

oxidation the crystals arc heated with a suitable amount 
(« 10 g) of YBCO buffer powder in a quartz tube con- 
nected to an external closed volume system. High purity 
oxygen gas (99.999 %) is supplied to the system and the 
pressure is controlled and monitored by use of high preci- 
sion absolute pressure gauges (MKS Baratron) with accu- 
racy and resolution better than 0.01 %. The crystals and 
the powder are wrapped in platinum foil to prevent reac- 
tion with the quartz tube. The closed volume system is 
made of ultra-high vacuum components and contained in 
a thermostatically controlled environment which allows 
for accurate determination of the oxygen pressure and 
oxygen uptake by the powder and the crystals. The wa- 
ter adsorbed in the system and the materials is removed 
prior to the preparation by use of a liquid nitrogen trap 
and heating the quartz tube to 300 °C. At this tempera- 
ture there is no reduction of the crystal and the oxygen 



The oxygen composition of the crystals was changed by 
use of a gas-volumetric equipmentElu. For reduction or 



equilibrium pressure is in any case sufficiently low that 
the oxygen does not condense in the trap. The desired 
oxygen composition x is usually established by pumping 
out or adding a suitable amount of oxygen gas at tem- 
peratures between 500 and 600 °C. During subsequent 
cooling the oxygen pressure is reduced to assure that x 
stays essentially constant. A final long time annealing 
may be performed at lower temperatures to obtain equi- 
librium between the powder and the single crystals and 
develop the oxygen ordering superstructures. For studies 
of the structural phase diagram a characteristic proce- 
dure to establish the superstructure is annealing at 80 
°C for 10 hours and cooling by 1 °C/hour to room tem- 
perature where the crystal is stored for more than one 
week before the measurements. 

If the starting oxygen composition of the buffer powder 
is known and it is assumed that the crystals are in equi- 
librium with the powder, the oxygen composition may 
be determined with an accuracy better than Ax = 0.02 
by use of the ideal gas law. The resulting oxygen com- 
position x has been compared with the known values of 
the oxygen equilibrium pressure determined by Schleger 
et alB, and full agreement has been established in all 
cases. Crystals prepared previously by the method have 
been examined by neutron diffraction technique and the 
oxygen composition x determined from crystallographic 
analysis of 375 unique reflections were found to be in 
full agreeiaaent with the values obtained from the gas- 
volumetrytD. 



B. Instrument 

The experiments were performed on a triple axis 
diffractometer at therhjgh-energy beam line BW5 at HA- 
SYLAB in Hamburger. The diffractometer operates in 
horizontal Laue scattering geometry and is equipped with 
a Huber 512 Eulerian cradle and a solid state Ge detec- 
tor. The insertion device is a high field wiggler with a 
critical energy of 26.5 keV at the minimum gap of 20 
mm. A 1.5 mm copper filter cuts the spectrum below 50 
keV, thereby minimizing the heat load on the monochro- 
mator crystal. The incident radiation with an energy in 
the range of 100 keV has a penetration depth of ~ 1 
mm in YBCO samples. For monochromator and ana- 
lyzer crystals either (2 0) SrTiOa crystals or (1 1 1) 
Si/TaSi2 crystalsEa were used. Both types of crystals 
had a mosaic spread of ~ 50" (arc seconds), resulting 
in a longitudinal resolution of 0.0075 A -1 at the (2 0) 
reflection of YBa2Cu306+a;. The transverse resolution is 
limited by the sample mosaicity, which was in the range 
of 0.05°-0.1° for our samples, corresponding to ~ 0.0015 
A^ 1 at the (2 0) reflection. The vertical resolution de- 
pends on the setting of the slits before and behind the 
sample. They were usually set to integrate the scatter- 
ing over a quarter of a reciprocal lattice unit, that is 0.40 
A^ 1 along the a and b axes and 0.13 A -1 along the c 



axis. The sample was wrapped in Al-foil and mounted 
in a small furnace, designed for use in an Eulerian cra- 
dle. The furnace temperature was stable within 1°C. An 
inert atmosphere of 0.3 bar Ar was introduced into the 
furnace to prevent oxidation of the crystals. From the 
gas-volumetric preparations it is established that the re- 
duction is negligible for temperatures below 300 °C, and 
we observed no changes in the structural properties that 
could be related to a change of oxygen composition dur- 
ing temperature cycling at temperatures below 250°C. 



C. Analysis of superstructure data 

The ortho-II superstructure reflections are well de- 
scribed by the scattering function 

5(q) = A/(l + (q h /T h f + {q k /T k f + (m/Ttff (1) 

where qi, i — h,k,l is the reduced momentum transfer 
and Ti the reduced inverse correlation length, related to 
the correlation length £j by £), = a/(27rl\) for the a di- 
rection and analogous along b and c. Equation |lj is a 
3D anisotropic Lorentzian raised to the power y. The 
scattering function S(q) has been convoluted with the 
resolution function, which is treated as a (5-function in 
the scattering plane and an integrating function in the 
perpendicular direction. The exponent y indicates the 
distribution of domain size, i.e. y = 1 points to an ex- 
ponential decrease of the pair correlations as for example 
in the ramified clusters typically for critical fluctuations 
above the transition temperature. The exponent y = 2 
may result from_a domain size distribution around an 
average value rjfj and, as mentioned in the Introduc- 
tion, the asymptotic behaviour for large q is in agree- 
ment with Porod's law for scattering from 3D finite size 
domains with sharp boundaries. Furthermore Bray has 
shown that the tail of the scattering function from a topo- 
logical defect of dimensiop-,r7i, in a system of dimension d 
is given by S(q) ex 2 }_ m E3. The relation between F and 

the peak width (HWHM=A) is A/^/v - 1. When full 
integration of the superstructure peak is performed per- 
pendicular to the scattering plane by relaxing the verti- 
cal aperture the in-plane scattering function derived from 
Eq. [His described by a Lorentzian to the power y' = y—\- 
The ortho-Ill, ortho-V and ortho-VIII superstructures 
are essentially 2D ordered, giving rise to significant over- 
lap of the peaks along KE3E3. In this case full integra- 
tion in the vertical direction cannot be obtained when the 
c axis is perpendicular to the scattering plane. For the 
2D ordered superstructures with finite domain size and 
sharp boundaries it is expected that the scattering func- 
tion in the ab plane should be a Lorentzian to the power 
y' = y = 3/2 because the integration along c is rather 
incomplete, whereas it should be a simple Lorentzian 
(y 1 = y — 1/2 = 1) when the c direction is in the scatter- 
ing plane and full integration along the a or b direction 
is performed. 



III. RESULTS 
A. The ortho-II superstructure 

The ortho-II phase is observed as a 3D ordered struc- 
ture with anisotropic correlation length at compositions 
0.35 < x < 0.62 and no other superstructure was found 
in this range of oxygen content. Especially at x = 0.35 
the q-space was surveyed unsuccessfully for correlations 
of the herringbone, \/2a x 2\/2a X c, or 2y/2a x 2\/2a x c 
type stmc turo s which-svere observed in other studies by 
electronl3't3E3, x-raye3 and neutron diffractionli-3. 

The largest correlation lengths are obtained for the 
ortho-II phase at x — 0.50. An example of the (2.5 0) 
superstructure reflection measured in crystal #1 is shown 
in Fig. [i] The crystal has been annealed by the ther- 
mal procedure marked 2 below. The room temperature 
properties of the oxygen superstructures are strongly de- 
pendent on the crystal quality and the thermal treat- 
ment. The four crystals labeled #1 - #4 in section [I A 
have been prepared with x — 0.50, annealed at 500 °C 
for 6 days, cooled by 10 °C/hour to 100 °C where they 
were annealed for 36 hours and then quenched to room 
temperature. The diffraction studies were performed 5- 
10 days later. With this thermal treatment we assume, 
on the basis of the results presented below and in Ref. 
ff4| , that all the crystals have reached the late state of 
ortho-II domain growth, and the influence of the differ- 
ent room temperature annealing times is considered to 
be small compared to the differences due to impurities 
and defects. 

Although it is possible to determine the impurity level 
in the crystal, it is often not known on which lattice sites 
the various impurities are located and what influence 
they have on the oxygen ordering. However, a crystal 
prepared as described in Ref. 



36 



using AI2O3 as crucible 
material resulted in a crystal with 6 mole % Al, which 

bund to be located on 
. This crystal did not 



by neutron diffraction studies wj 
the Cu(l) site in the basal plan- 
show ortho-II superstructure ordering when it was pre- 
pared with x = 0.50. The nature and influence of lattice 
defects are equally difficult to quantify. One way to mea- 
sure the overall quality of the crystal lattice is the mo- 
saicity, i.e. the width of the rocking scan of the sample. 
Fig. g shows the HWHM of h, k and I scans obtained 
at room temperature of the ortho-II superstructure re- 
flection versus the mosaicity for the four crystals #1 - 
#4. The h and I scans are measured in the ac scatter- 
ing plane, and the scan along k in the ab plane. In both 
cases a full integration of the respective vertical widths, 
Afc and Al is performed. It is immediately obvious that 
crystal ^3 prepared with low purity chemicals have the 
largest mosaicity and the broadest ortho-II peaks. Thus, 
the purity of the chemicals is crucial for the development 
of large ortho-II domains. For crystals #1, #2 and #4 a 
linear relation between mosaicity and the HWHM of the 
ortho-II superlattice peak is found, with small deviations 



of the width along h. With the mosaicity as a criteria of 
crystal quality crystal #1 is the most perfect one. This 
is corroborated by magneto-optic studies of the magnetic 
flux flow in the crystals. Crystal #1 is the only one that 
shows flux flow instability, which is considered to be a 
signature of very high crystal qualitycj. 

The influence of the thermal treatment has been stud- 
ied in crystal #1. After the preparation for x — 0.50 
mentioned above the crystal has been annealed the fol- 
lowing six ways: 

1. 70 days annealing at 80° C, cooling to room tem- 
perature in steps of f °C/hour 

2. 5 hours annealing at f 00° C, quenched to room tem- 
perature and stored for 10 days 

3. quenched from 170°C to room temperature and 
stored at room temperature for 97 days 

4. cooled down from 170°C to room temperature in 
steps of 10°C every 10 minutes 

5. cooled down from 170°C to room temperature with 
4°C/minutes 

6. quenched from 170°C to room temperature within 
3 minutes 

In Fig. the normalized peak intensity is plotted versus 
the inverse correlation lengths measured at room tem- 
perature at the (2.5 5) superstructure reflection for 
the differently treated samples. The normalization of 
the peak intensity is made relative to the background. 
As mentioned above, the h and I scans are measured in 
the ac plane, the scans along k in the ab plane, and a full 
integration along the vertical widths is performed. The 
relation between the measured peak intensity and the 
peak widths in all three directions follows a quadratic 
dependence as indicated by the lines. This shows that 
the ratios between the line- widths Th/^k and Ti/Tk are 
independent of the thermal treatment. Further, since 
the measured peak intensity, I°^ k , includes an integral 
over the direction perpendicular to the scattering plane 
we establish that also the total integrated intensity, Ii„t , 
is independent of thermal treatment. This is easily seen 
from Fig. and the following relations, where the inte- 
gration is assumed to be along k: 

*int &- -LpeakL /jl fc-1- / 
Ipeak ^ Ipeak^k OC Iint/{^hTl) OC lint /T \ 

Accordingly, only the correlation length of the super- 
structure and thereby the characteristic domain size de- 
pend on the sample treatment. This indicates that the 
finite size domains have internal thermodynamic order 
and fill the crystal. Studies of the time dependent oxygen 
ordering following a temperature quench from the ortho-I 
into the ortho-II phase at this composition confirm that 



the integrated intensity depends only on the ptempera- 
ture. These results will be published elsewhereci 

It is instructive to consider the temperature variation 
of the ortho-II structure for crystal #1 prepared by ther- 
mal treatment 1. The phase transition from the ortho-II 
phase into the ortho-I phase was studied by means of 
w-scans at the (2.5 5) reflection during initial heating 
and at subsequent cooling within one hour. The result is 
shown in Fig. 0. Clearly, the peak intensity is lower and 
the u> width is larger on cooling than on heating. How- 
ever, the integrated intensity calculated as I°^, s ak x to 2 is 
found to be the same during heating and cooling. Since 
the results of the ortho-II superstructure ordering indi- 
cate that internal superstructure order is established in- 
side the finite size domains it is appropriate to define a 
transition temperature, Ton- Several criteria may be 
used. Firstly the variation of the peak intensity of the 
superstructure reflection plotted in the top part of Fig. 
[I shows an inflection point at 95° C, as determined by 
the minimum of the normalized slope (N.S.) of the peak 
intensity, plotted in the inset. The inflection point of the 
peak intensity indicates the cross-over from static order 
to critical fluctuations, i.e. the transition temperature 
Ton- Secondly the onset of the line broadening of the 
superstructure reflection also marks the transition tem- 
perature. The temperature dependence of the line width 
above the transition can be well described by the critical 
exponent v — 0.63 for the 3D Ising model: 



A(T) = A±|T-T C | 



(2) 



where Aq are the slopes below and above the transi- 
tion temperature. This behavior is shown in the middle 
part of Fig. H. From the fit to the heating data in the 
critical region we find Ton = 95 °C. Thirdly the line 
shape changes from approximately a Lorentzian squared, 
i.e. y' = §, at room temperature to a simple Lorentzian 
at the transition temperature: Ton = 95°C shown in 
the bottom part of Fig. || A line shape described by 
a simple Lorentzian is characteristic for critical fluctua- 
tions above the transition temperature. During cooling a 
drastic slowing down of the ordering process is observed 
at temperatures close to the transition temperature as 
seen in the variation of the peak intensity, which starts 
to deviate from the heating data at ~ 105°C. 

The variation of the peak intensity and the peak width 
with temperature leads to the distinction of three areas 
during the heating cycle. Between room temperature and 
50°C both the peak intensity and the width of the super- 
structure reflections are constant, i.e. both the domain 
size and the integrated intensity, which is a measure of 
the order parameter of the ortho-II phase, are constant 
within the time period studied. Between 50°C and 95°C 
the width of the superstructure reflections is still con- 
stant, but the peak intensity decreases with increasing 
temperature, indicating that the ortho-II order inside 
the domains and thereby the number of oxygen atoms 
ordered in alternating full and empty chains start to de- 
crease. Finally, the increasing width and the decreasing 



intensity in the temperature range above 95°C indicate 
the area of critical fluctuations above the transition tem- 
perature. This range is well described by the critical ex- 
ponents. In contrast, the behavior below the transition 
temperature shows substantial deviations from what is 
expected from a regular second order phase transition, 
where a long range ordered phase should be formed. 

The investigation of the temperature dependence of 
the ortho-II phase at x = 0.42 exhibits exactly the same 
behavior as found at x = 0.50 with a small shift in the 
transition temperature. Thus a stoichiometric oxygen 
content is only of minor importance for the behavior at 
the ortho-I/ortho-II phase transition. However, the peak 
widths, Ah = 0.031(1), Afc = 0.010(1), Al = 0.14(1), are 
significantly larger than found for the high purity crystal 
#1 with x = 0.50, cf. Figs. | and §. 



B. The ortho-Ill structure 

The ortho-Ill phase is found at the oxygen composi- 
tions of x — 0.72, 0.77 and 0.82. A crystal prepared with 
x = 0.87 showed no sign of any oxygen ordering. At 
these oxygen compositions the ortho-Ill phase is formed 
by the sequence (110) of two full chains and one empty 
chain. Accordingly, the size of the unit cell is tripled 
along a and the diffraction pattern shows two superstruc- 
ture reflections along h between the fundamental Bragg 
peaks. As shown in Fig. the ortho-Ill super lattice 
peaks are well defined in the ab plane, but like all super- 
structures due to oxygen ordering in YBCO, broadened 
due to finite domain sizes. In contrast to the ordering 
in the ab plane, the /-dependence of the diffracted inten- 
sity shows only a broad modulation, with a peak width 
corresponding to more than one reciprocal lattice unit. 
This /-modulation is characteristic of the ortho-Ill struc- 
ture and has bees-, found in all samples exhibiting the 
ortho-Ill phaseE3'CJ'E!3, which indicates that the order- 
ing of oxygen atoms takes place in the a6-plane, whereas 
different planes are only weakly correlated. Thus, in con- 
trast to the ortho-II structure, which is 3D ordered, the 
ortho-Ill phase is essentially a 2D ordered superstruc- 
ture. As mentioned in Section II C a simple Lorcntzian 
is expected for scattering from finite size domains in a 2D 
system when the integration of the peaks is performed ei- 
ther along the a or the b direction (y' — y — 1/2 = 1), 
whereas a Lorentzian to the power y' = y = 3/2 is ex- 
pected for integration along I. However, the peak shapes 
along h and fc are well described by Lorentzians, inde- 
pendent of which component is integrated in the vertical 
direction. Attempts to include a variable power y did not 
improve the fits when the I direction was vertical. The 
smallest widths,, which have been reported previously by 
Schleger et al!E3 are found in the x = 0.77 crystal with 
Ah = 0.031(1) and Afc = 0.0090(2). 

One example of the temperature dependence of the 
ortho-Ill phase is shown in Fig. @ for the crystal with 



composition x — 0.72. Here the (8/3 5) reflection was 
scanned along h at various temperatures. Similar to the 
transition of the ortho-II phase the peak intensity and 
peak width are frozen at temperatures smaller than 35° C. 
Above this temperature critical fluctuations are observed. 
Fitting the temperature dependence of the peak width 
a critical exponent of v = 0.92(8) is obtained, with a 
transition temperature of Tom = 48(5)°C. This value 
for the critical exponent is in good agreement with the 
theoretical value of v = 1 for the 2D Ising model and 
confirms the 2D character of the ordering. 



C. Ortho-V 

The investigation of a crystal prepared with the oxy- 
gen composition of X = 0.62 shows a mixture of ortho- 
II and ortho-V phase at room temperature. This is re- 
vealed by the observation of diffuse peaks at positions 
of h = 2.4, 2.5 and 2.6 as shown in Fig. m. The peak 
at ft. = 2.5 results from the ortho-II structure, and the 
peaks at h = 2.4 and h = 2.6 are consistent with a unit 
cell which is enlarged five times in the a direction, i.e. 
the ortho-V structure. The two small peaks seen in the 
/i-scan in Fig. [?] at h—2.23 and /i=2.83 are an Al-powder 
line and possibly a grain of an unknown phase oriented 
with the lattice, respectively. The hump at /i=2.83 has 
also been observed when the same crystal was prepared 
with other oxygen stoichiometrics (compare with Fig. O 
and Ref. |26| ). A similar diffraction pattern, consistent 
with a mixture of ortho-II and ortho-V has been ob- 
served in all (h I) scans performed with 1 < h < 4 and 
Z=0,3,5,6,7 (8 scans in total). However, none of these 
scans showed a peak at position Q = (y 1). This is 
explained by the structure factor calculations of the su- 
perlattice peaks from the ideal ortho-V ordering sequence 
(10110) shown in Ref. 47. The intensities of the peaks 
at (1/5 0) are indeed much smaller than the ones at 
(2/5 0) and (3/5 0). However, this model takes 
into account only the oxygen order, and, as discussed in 
Section |, the superlattice peaks are caused by both the 
oxygen order and the cation displacements. These dis- 
placements and the pronounced disorder may change the 
intensities and reduce them further. 

Due to the heavy overlap of the peaks from the two 
phases it is difficult to determine the peak shape and 
width. However, analysis of the ortho-II and ortho- 
V peaks using Lorentzian profiles gave the following 
HWHM in reciprocal lattice units at room temperature: 
Ah = 0.040(28), Afc = 0.0078(16), A/ = 0.12(3) for 
ortho-II, and Ah = 0.058(10), Afc = 0.0096(19)for ortho- 
V. The scan along I at (2.5 I), shown in Fig. [?], exhibits 
the intensity modulation well known for 3D ordering in 
the pure ortho-II phase, but the heavy overlap of the 
peaks prevents that the 2D short range type of mod- 
ulation expected for the ortho-V peaks along I can be 
determined independently. 



The temperature dependence of this mixed phase was 
measured by ft-scans between the (2 0) and the (3 0) 
Bragg reflections and the diffraction pattern was fitted to 
three Lorentzians with fixed positions at 2.4, 2.5 and 2.6. 
Looking at the measurement of the phase transition of 
this mixed phase, shown in Fig. || one observes that the 
ortho-V correlations disappear between 50° C and 70° C 
and at the same time the ortho-II gains intensity. Also 
during the cooling cycle the ortho-II correlations dom- 
inate the diffraction pattern. The ortho-II correlations 
disappear at approximately 110 °C. Both factSptpffcther 
with our knowledge about the ordering kineticsauJ indi- 
cate that during the cooling process the ortho-II phase is 
stabilized at higher temperature than ortho-V and with 
rather fast ordering kinetics. Then at lower temperature 
the ortho-V phase becomes stable, but due to the slow 
ordering kinetics at this lower temperature the ortho-V 
domains do not form within the one hour time period of 
the experiment. 



D. Ortho-VIII 

Figure |9j shows h, k and I scans for the oxygen com- 
position x — 0.67. The left part with h scans along 
(h 0) with 2 < h < 3 reveals diffuse superlattice peaks 
at h = 2.382(4) and h = 2.627(3). The peak positions 
and profiles have been fitted to two Lorentzians giving 
a HWHM of Ah = 0.053. The middle part shows that 
the peaks are also localized in the transverse direction 
with a width of Afc = 0.013(2). The modulation of the 
intensity for a scan along I (right part of the figure) has 
a similar q dependence as the corresponding scan for the 
ortho-Ill phase, shown in Fig. ||. Thus, there are no 
well-defined peaks along I, indicating essentially 2D or- 
dering with substantial disorder in the stacking of full 
and empty chains along the c direction. Similar superlat- 
tice peaks have been observed at positions in reciprocal 
space of (h 3) and (h 5) with 2 < h < 3. The su- 
perstructure peak positions at modulation vectors with 
nh m = 0.382 and 0.627 are close to the expected values 
nh m = 3/8 and 5/8 for a superlattice with a unit cell 
of 8a x 6 x c, i.e. the ortho-VIII phase. The expected 
sequence of full and empty chains of the ideal ortho-VIII 
structure is (11010110). Calculating the intensities of the 
superlattice peaks for this ideal case one finds that the 
observed peaks at nh m = 3/8 and 5/8 are the strongest, 
the peaks at nh m = 2/8, 4/8 and 6/8 are about one order 
of magnitude smaller, and the ones at nh m — 1/8 and 
7/8 are about two orders of magnitude smaller (compare 
with the presentation in Ref. E^). Due to the weak order- 
ing it is unlikely that the smaller superlattice reflections 
can be observed. 

The temperature dependence was observed from h- 
scans at the (2j^ 3) peak position and the results are 
shown in Fig. [fc]. The onset of broadening of the super- 
structure peaks takes place at Tqviii = 42(5)°C. The 



temperature dependence of the peak width above the 
transition temperature is described by the critical expo- 
nent of v = 0.79(3) as shown is the middle part of Fig. 
[Uj . This value is between the exponent of 0.63 for the 
2D and 1 for the 3D Ising model. Another interesting 
feature of this phase transition is revealed by the inspec- 
tion of the peak position. When the temperature ex- 
ceeds 50 °C the peak position changes continuously from 
h = 2.372 to h = 2.4 which corresponds to the position 
of the ortho-V phase. Above 90°C the peak shifts grad- 
ually to ft. = 2.33 at 150°C, the location of the peaks of 
the ortho-Ill structure. Upon cooling the data are repro- 
duced down to 75° C, at lower temperatures the intensity 
is significantly reduced and the structure freezes into the 
ortho-V phase. 



E. The oxygen ordering phase diagram 

From the transition tempejjstijrcs obtained in the 
present and previous studies&EBu using hard x-ray 
diffraction and the same type of crystals we may estab- 
lish phase lines for the oxygen superstructure ordering. 
Combining these data with the transitions temperatures, 
To i, of the phase transition from the tetragonal to the 
orthorhombic ortho-I phase, obtained by neutron pow- 
der diffraction!!!, we may construct the structural phase 
diagram of oxygen ordering in YBCO, shown in Fig. 
O . In the figure is also included the phase transition 
temperatur.es, Toi and Ton predicted by Mpiate Carlo 
simulationseS based on thepASYNNNI modelij with ab 
initio interaction parmetersEj. 

The only true equilibrium structures are the ortho- 
I phase and the tetragonal phase, all superstructures 
formed by oxygen ordering do not show long range order- 
ing. Within the temperature range studied the tetrago- 
nal phase is the only one observed for x < 0.35. Be- 
low the tetragonal to orthorhombic phase transition tem- 
perature the 3D ordered ortho-I phase always develops, 
and it is the only structure observed for x > 0.82. For 
0.35 < x < 0.62 the 3D short range ordered ortho-II 
phase is the only stable superstructure. Similarly, a sin- 
gle phase ortho-Ill structure with 2D finite size ordering 
is observed for 0.72 < x < 0.82. At intermediate compo- 
sitions a mixed phase of ortho-II and ortho-V is found at 
x = 0.62, and ortho-VIII is found at x — 0.67 in crystals 
that have been slowly cooled to room temperature as de- 
scribed in Sec. II A . Both the ortho-V and the ortho-VIII 



structures are essentially 2D ordered and have finite size 
ordering. During heating the ortho-V structure trans- 
forms into ortho-II and it does not recover on cooling 
within one hour. Above room temperature the ortho- 
VIII structure transforms gradually first into ortho-V and 
then into ortho-Ill. On subsequent cooling the ortho-V 
superstructure is recovered and remains stable within the 
one hour time period of the measurements. 

The line shape of the superstructure reflections is in 



most cases well described by a simple Lorentzian (y = 1). 
Only for the ortho-II phase between 0.42 < x < 0.62 a 
Lorentzian squared shape (y = 2) is found. At the low 
oxygen side of the ortho-II phase x < 0.36 the small 
peak to background ratio (see bottom part of Fig. |l2] ) 
does not permit the determination of the exponent of the 
Lorentzian. The domain size of the superstructures de- 
pends strongly on the crystal quality and the annealing 
times. However, for high quality crystals that have been 
annealed by the standard proced ure fo r studies of the 
phase diagram (described in Sec. II A) we expect that 



the domains are at the late state of growth and there- 
fore only weakly time dependent (c/. Fig. 0, thermal 
preparations 1 and 2, and Refs. || and p4|). On this basis 



we consider the results presented in Sec. |I| of the peak 
widths measured at room temperature after the initial 
thermal preparation as saturation values. The HWHM 
of the superlattice peaks measured along the three axis 
of reciprocal space as function of oxygen composition is 
depicted in Fig. |l2| (top). The parallel lines (guides to 
the eye) in the logarithmic plot observed in the ortho-II 
phase as well as in the ortho-Ill phase show that the ra- 
tio of the anisotropy is constant within a given structural 
phase. For the ortho-II phase we find the following ratios 
of the inverse correlation lengths at room temperature: 
T h /T k = 2.7(6) and r;/r fc = 15(2). The ab plane ratio 
seems to be independent also of the type of structure, 
since the ratio for the ortho-Ill phase: Th/Tk = 2.9(4) is 
in good agreement with the value of the ortho-II phase. 
This implies that the domain pattern in the ab plane 
scales in both the ortho-II and the ortho-Ill phases, and 
in ortho-II the scaling is extended to 3D. The peak in- 
tensities cannot be compared directly because different 
crystals and instrumental settings have been used. How- 
ever, the peak intensity normalized to the background, 
shown in the bottom part of Fig. [12[ is essentially an 
independent parameter of the ordering properties. From 
this normalized peak intensity and the HWHM data it is 
clear that the optimal superstructure order parameter is 
found close to x = 0.55. 

The oxygen composition x of all the ordered phases de- 
viates systematically from the ideal composition of these 
phases. For example the longest correlation lengths for 
the ortho-II phase are likely to be at x ~ 0.55. Unfortu- 
nately no data points are available at this composition. 
Theoretically one would expect the best ortho-II order 
for x = 0.50. This deviation is even more significant for 
the ortho-Ill phase, which is expected at x = 0.67, but 
observed around x ~ 0.77. Thus, the deviation from the 
ideal composition increases with increasing chain density 
and the amount of oxygen atoms occupying sites on the 
empty chains at room temperature can be estimated to 
be about 10 % for ortho-II and 30 % for ortho-Ill. 



IV. DISCUSSION 



A. Experimental results 

It has been known for several years that the ortho-II 
and ortho-Ill superstructures are bulk structural phases 
of finite size domains. Several other superstructures have 
been suggested, mainly from electron microscopy. In the 
present paper we have shown that also the ortho-V and 
ortho-VIII correlations observed by electron microscopy 
result from bulk structural ordering, but we found no 
evidence for the ortho-IV phase. However, we recognize 
in particular the early .-electron microscopy results ob- 
tained by Beyers et aZ.L2l which are in close agreement 
with our room temperature data. Beyers et al. observe 
the ortho-II and ortho-Ill superstructures in the same 
composition range as in our studies. Furthermore, they 
found co-existence of ortho-II and ortho-V at x = 0.65, 
and a structure similar to the ortho-VIII phase, which 
they call a '(0.37 0)' structure, at x = 0.71. 

Beyers et al. attributed the clear disagreement be- 
tween the observed oxygen compositions and the stoi- 
chiometrics of the ideal superstructure phases to gradi- 
ents in the oxygen content of the sample, which might be 
different on the surface and in the bulk material. In our 
experiment such differences can be ruled out. We con- 
clude that this deviation is an intrinsic property of the 
oxygen ordering mechanism. It is possible that the phase 
lines between the superstructure phases are in fact tilted, 
and only at zero temperature the ideal oxygen stoichiom- 
etry of the superstructure phases is found. However, this 
will never happen because the oxygen ordering kinetics is 
very slow at the temperatures where the superstructures 
become stable, and the movement of Cu-0 chains freeze 
effectively below approximately 40 °C. 

Beyers et al. interpret the mixing of ortho-II and 
ortho-V phases at x = 0.65-as a phase separation, which 
leads to the 60 K plateauI3. Our investigation of the 
temperature dependence together with the studies of the 
ordering kineticsuO may lead to a different conclusion. 
During the cooling of a sample with an oxygen content 
of x — 0.62 (in the case of Beyers et al. x = 0.65) oxy- 
gen starts to order in the ortho-II phase. The relatively 
high temperature enables a fast growth of ortho-II do- 
mains. At lower temperature the ortho-V phase becomes 
stable, but now at temperatures just above room temper- 
ature where the growth of ortho-V domains is slow. A 
full transformation into the ortho-V phase cannot be pre- 
cluded but it is very time consuming. Therefore, we sug- 
gest that domains of the complex ortho-V superstructure 
start to grow inside the ortho-II structure and a mixed 
phase, rather than phase separation, results. 

From studies of the oxygen ordering properties it has 
become clear that the finite size of the ortho-II super- 
structure results from formation of anti-phase bound- 
aries that limit the domain growth due to slow kinetics 
of moving long Cu-0 chains. The reason for this has 
been discussed by Schleger et alB, and it was speculated 
that random fields introduced by impurity defects in the 



crystal stabilize the anti-phase domain walls and prevent 
formation of long range order. This is corroborated by 
the present.-studies and additional studies of the order- 
ing kineticalj. However, the observation of superstruc- 
tures extending over eight unit cells shows the impor- 
tance of long range interactions for the ordering mech- 
anism. That these long range interactions play a sig- 
nificant role for the finite size ordering has recently been 
established by model simulations!^, and will be discussed 
further below. For the ortho-Ill, ortho-V and ortho-VIII 
superstructures the small 2D domains indicate that the 
ordering resembles a random faulting sequence of ortho- 
II and ortho-Ill. Khachaturyan and MorrisEi have sug- 
gested that this is a likely ordering scheme, and they 
have calculated structure factors which are qualitatively 
similar to those observed at room temperature in our 
experiments. However, the fact that the ortho-V and 
ortho-VIII superstructures only appear when they are 
slowly cooled indicates that the long range interactions 
tending to form these superstructures become effective 
at low temperatures, but the slow oxygen ordering ki- 
netics for movement of long Cu-0 chains prevent that 
well-defined domains are formed. As mentioned in Sec. 



II C we would expect a diffraction profile of a Lorentzian 
to the power y' = y = 3/2 from 2D domains with sharp 
boundaries when the integration along the c axis is in- 
complete. The observation that all the superstructure 
peaks of the ortho-V, ortho-VIII and ortho-Ill peaks are 
described by Lorentzian profiles suggests that these su- 
perstructures have a more fuzzy type of boundaries than 
the ortho-II domains. 

Generally, there is significant hysteresis in the su- 
perstructure ordering when the temperature is cycled 
through the phase transitions. The ortho-II and ortho- 
Ill superstructures are re-established during cooling from 
the ortho-I phase within one hour. However, the ortho- 
V phase (mixed with ortho-II) and the ortho-VIII do 
not recover during cooling within this short time pe- 
riod. Instead, the less complex superstructures, ortho- 
II and ortho-V develop, respectively. It is obvious that 
the superstructure ordering does not represent equilib- 
rium phases, and it cannot be precluded that more com- 
plex superstructures may be formed by very long time 
annealing at an appropriate temperature or in crystals 
that are even more perfect than the present ones. Ac- 
cording to Ostwald's step rule for phase transformations, 
metastable phases may be formed, before the system fi- 
nally transforms into the stable phase, as long as nucle- 
ation centers with a similar structure like the metastable 
phases are present. In our case the ortho-II and ortho-Ill 
phase might be nucleation centers for the ortho-V phase, 
which in turn, at X = 0.67, is metastable and serves as a 
nucleation center for the ortho-VIII phase (see figure 10). 
Thus, although we have been able to define unique transi- 
tion temperatures, at least for the ortho-II superstructure 
phase, it is questionable whether we have established a 
phase diagram in the usual sense. This may explain why 
the phase diagram does not comply with Gibb's phase 



rule. 

The 3D ordered superstructures with unit cells 2 V2 a x 
2\[2a x c and \/2a x 2\[2a x c, the so called her- 
ringbone structure, have been observed by electron mi- 
croscopy, and one groupJias reported on these structures 
by neutrontj and x-rayc3 diffraction techniques on sin- 
gle crystals with composition x = 0.35. However, no 
other experiments with bulk structural-techniques, could 
confirm these results. Bertinotti et alE3 and Yakhou et 
alx3 have shown that the reflections of the herringbone 
type can be assigned to BaCu304 grains in the crystals. 
Krekels et alBiE^ attribute the 2\[2a x 2\[2a x c struc- 
ture to distortions of tha-CuOs pyramids in the Cu02 
planes, and Werder et alEB suggest that they could re- 
sult from ordering of copper and barium vacancies in the 
lattice. The consensus from these and several other stud- 
ies is therefore that the 2\[2a x 2\J~2a x c and the her- 
ringbone type structures are not oxygen ordering super- 
structures in YBCO. If they were, it is peculiar that they 
have 3D long range order while the Cu-0 chain ordering 
develops only finite size domains. Also, we have found 
no evidence of them at any composition x in the present 
hard x-ray diffraction studies on carefully prepared high 
quality single crystals. 



B. Significance for superconductivity 

The significance of the oxygen ordering for charge 
transfer and superconductivity is obvious from many 
studies. Chemical bond considerations^pqmbined with 
structuralEZI and spectroscopic studiesEfEa have shown 
that the basal plane copper in undoped YBa2Cu3 06+£c, 
(x = 0), is monovalent and that simple oxygen 
monomers, i.e. Cu-O-Cu, will not give rise to charge 
transfer. However, charge transfer is observed foiuLarger 
x where Cu-0 .chains are formed. Cava et alEU and 
Tolentino et alE3 have established that an increasing 
amount of oxygen give rise to a charge transfer to the 
Cu02 planes that is in good agreement with the well- 
known plateau variation of T c with T c = 58 K around 
x = 0.5 and T c = 93 K close to x = 1. Relating the oxy- 
gen ordering to the variation of T c observed e.g. by Cava 
et alrli we find that the 58 K plateau is identical to the 
stability range of the ortho-II superstructure, the rise in 
T c from the 58 K to the 93 K plateau takes place at val- 
ues of x where the ortho-V/II, ortho-VIII and ortho-Ill 
structures are found, and the 93 K plateau coincide with 
the oxygen compositions of the ortho-I phase. 

The significance of the ortho-II ordering for supercon- 
ductivity has .-been shown directly by Veal et ala and 
Madsen et alii. Both groups have shown that T c is sig- 
nificantly reduced just after a quench and increases with 
time towards the equilibrium values with a thermal acti- 
vated time constant. The conclusion, that may be drawn 
from these experiments and the present structural data, 
is that the formation of ortho-II superstructure is deci- 
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sive for the charge transfer and T c . When the sample is 
quenched from temperatures above Ton, and even from 
the tetragonal phase, it is only the time used to quench 
it into the ortho-II phase that matters. A time depen- 
dent increase of T c is observed at annealing temperatures 
down to 250 K. Most likely this temperature is the lower 
limit for local oxygen jumps which dominates the oxygen 
ordering at the very early time. It is unlikely that the 
domain wall separating anti-phase domains are mobile at 
250 K. 



C. Relation to model calculations 

There has been many theoretical studies of the oxy- 
gen ordering in YBCO and attempts to correlate the 
structural ordering with the electronic properties and su- 
perconductivity. These include phenomenoloffiicaL rela- 
tions between ordered oxygen domains and T c E2tE!j, and 
electron band structures calculated from oxygen chain 
configurations estimated ad hocEi or derived from model 
studiesEJ. More realistic and elaborate models, where 
the electronic degrees of freedom from the strongly cor- 
related electron system has been included in combina- 
tion with the-oxygen ordering properties, have also been 
consideredE-3~E3. The aim is clearly to understand details 
of the local oxygen ordering properties which are impor- 
tant for the electronic structure and the charge transfer, 
but difficult to obtain directly from experiments. The 
predictive power of these model studies for the structural 
and electronic properties is strongly related to their abil- 
ity to reproduce the experimental findings, as presented 
in the present structural studies. 

Most of the structural models are based on local ef- 
fective oxysaj-oxygen interactions in a 2D lattice gas 
formulationca'EZl. These models do not take into ac- 
count long range interactions like strain effects and will 
therefore not reproduce the twin domain formation ob- 
served experimentally at the onset of the ortho-I order- 
ing. Models including such long range interactions have 
been considered^. It has also been suggested that the 
diffuse scattering results from generation of (1 0) in- 
terstitial plane defects [that order by forming a Magneli 
type homologous seriescj. However, the oxygen ordering 
observed experimentally has a predominant 2D character 
related to the CuO^ basal plane, that is (0 1) planes. It 
takes place inside the ortho-I twin domains and has do- 
main sizes that are usually much smaller than the twin 
domains. Strain effects and long range Coulomb type in- 
teractions are therefore of little significance for the oxy- 
gen superstructure formation but they do play a role for 
the mesoscopic ordering properties. 

The simplest model, that accounts for many elements 
of the oxygen ordering properties, like the formation of 
Cu-0 chains and the presence of the tetragonal, ortho-I 
and ortho-II phases, is the so-called ASYNNNL (Asym- 
metric Next Nearest Neighbor Ising) modelEJ. The 



ASYNNNI model is a 2D lattice gas (or Ising) model 
with effective oxygen-oxygen pair interactions, that are 
assumed to be independent of temperature and compo- 
sition x. The interactions parameters include a strong 
Coulomb repulsion V± between the oxygen atoms on near- 
est neighbor sites and an attractive covalent interaction 
V<i between oxygen atoms that are bridged by a Cu atom. 
These two interactions locate the oxygen on the b axis 
(the 0(1) site) and prevent significant oxygen occupa- 
tion on the a axis (the 0(5) site) in the orthorhom- 
bic phases at moderate temperatures and compositions 
x > 0.35. A weaker effective repulsive Coulomb type in- 
teraction V3 between oxygen atoms that are next nearest 
neighbors and not bridged by a Cu atom stabilizes the 
ortho-II superstructure. The ASYNNNI model accounts 
quantitativelyca for the temperature and composition de- 
pendence of the experimental structural phase,-transition 
between the tetragonal and the ortho-I phasesfcl (see Fig. 
|ll| ) by use of interaction parameters, which ace consistent 
with values obtained by Sterne and WilleEJ from first 
principles total energy calculations: Vi/ks = 4278 K, 
V 2 /k B = 1488 K and V 3 /k B = 682 K. It also predicts the 
existence of the ortho-II phase, but it cannot account for 
the additional superstructure phases, ortho-Ill, ortho-V 
and ortho-VIII. Moreover, it predicts long range order of 
the superstructure phases, which has never been obtained 
experimentally, and it cannot account quantitatively for 
the ortho-I to ortho-II phase transition temperature. 

Extensions of the ASYNNNI model have been sug- 
gested to account for the shortcomings. These include an 
effective 3D interaction with a nearest neighbpi; attrac- 
tive interaction along c that is V4 ~ — 0.02ViEia, effects 
of electronic degrees of freedom in the Cu-0 chain struc- 
ture, as mentioned aboveEji£3, and 2D C oulo mb type 



interactions of longer range than V2 and Vj 



For the 



2D ordering it has been argued that a single additional 
interaction parameter for oxygen atoms that &~£&,2a apart 
and not bridged by copper should be sufficients. This is 
corroborated by an estimate of screened Coulomb poten- 
tials which shows that the interaction between oxygen 
atoms separated by 2a is of the order of V5 = 0.021^, 
and it decays rapidly for larger distances. At the tem- 
perature where V5 becomes effective Cu-0 chains have 
already been formecS and it will act as an effective inter- 
action between chains rather than between oxygen pairs. 
The V5 interaction stabilizes the ortho-Ill phase by con- 
struction but it is not expected to account for the ortho- 
V and ortho-VIII phases. The influence of including ef- 
fective Coulomb type Cu-0 chain interactions extending 
beyond a and 2a has been studied analytically in the 
framework of a ID Ising modele! Here a sequence of 
branching phases develops for T — > in order to comply 
with the Nernst principle and stoichiometric phases at 
different compositions x. However, for the YBCO system 
it is expected that the interactions of range beyond 2a 
play a role only at low temperatures where effectively the 
structural ordering is frozen. Also, the projection to a ID 
system requires that rather long Cu-0 chains are formed, 
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and recent Monte Carlo simulations have shown that fi- 
nite chain lengths may result when the ASYNNNIjnodel 
is extended by the V5 interaction, even for T — v Cc3. On 
the other hand, the ASYNNNI model extended this way 
predicts that the Vg parameter is sufficient to establish 
not only short range correlations of ortho-V and ortho- 
VIII but also the ortho-II and ortho-Ill superstructures 
do not develop long range order, as observed experimen- 
tally. The finite size ordering of the superstructures is 
therefore not necessarily a consequence of impurities or 
defects that pin the domain walls, but may be an in- 
trinsic disordering property. Experiments on even more 
perfect crystals than used in the present study could sup- 
ply additional information about the influence of defects 
on the ordering properties. In further agreement with the 
experiments, the ASYNNNI model extended with the V5 
as well as the 3D V4 parameters predicts a significant 
suppression of the Ton ordering temperature relative to 
the Toi temperature which the original version failed to 
do (see Fig. |ll]). It is therefore, a promising model for 
analysis of experimental results and credible predictions 
about the local oxygen ordering properties. So far the 
theoretical phase diagram including the V4 and V5 inter- 
actions has not been determined. The present experi- 
mental results can be used as a guide to further model 
studies. Here it is interesting to note that our data show 
that the ratios of the correlation lengths are essentially 
independent of the oxygen stoichiometry in the ortho- 
II and the ortho-Ill phases. A comparison between_this 
result and mean field predictions of the peak widthsc3 in- 
dicates that the ASYNNNI model interaction parameters 
are independent of x. This assumption has been a major 
objection against the validity of the ASYNNNI model. 



nal thermodynamic equilibrium are formed. The domain 
size observed on cooling from the ortho-I phase within 
one hour is significantly reduced compared to the value 
obtained by long time annealing. The observation of 
ortho-V mixed with ortho-II, and ortho-VIII superstruc- 
tures shows that these superstructures are bulk proper- 
ties, and that Coulomb interactions beyond next-nearest- 
neighbors become effective close to room temperature. 
The ordering of the ortho-V and ortho-VIII superstruc- 
tures does not reproduce when the sample is cooled from 
the ortho-I phase within one hour, and it can not be 
precluded that additional superstructure phases may be 
formed by careful annealing of high quality single crys- 
tals. Therefore, although an unambiguous criterion has 
been identified for the ordering temperatures of the finite 
size ortho-II and ortho-Ill superstructures, the resulting 
'phase diagram' is not an equilibrium phase diagram in 
the usual sense. 
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V. CONCLUDING SUMMARY 



High energy X-ray diffraction has proven to be a 
unique tool for studies of oxygen ordering properties in 
the orthorhombic phase of YBCO. Chain ordered super- 
structures of the ortho-II, ortho-Ill, ortho-V and ortho- 
VIII types have been observed in high quality single crys- 
tals with this bulk sensitive technique. None of the super- 
structures develops long range order. Only the ortho-II 
phase is a 3D ordered superstructure with anisotropic 
correlation lengths. The ortho-II correlation lengths ob- 
served at room temperature depend on the oxygen com- 
position (optimal for x — 0.55), crystal perfection and 
thermal annealing. All other superstructures have 2D 
character with ordering only in the a6-plane. The ra- 
tio of the ab plane correlation lengths is essentially in- 
dependent of the oxygen composition and whether the 
ordering is ortho-II or ortho-Ill. The transition temper- 
atures of the superstructures are between room temper- 
ature and 125°C. The ordering properties resulting from 
thermal cycling through the Ton and the Tom ordering 
temperatures show that finite size domains with inter- 
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FIG. 1. (2.5 0) superstructure reflection of the ortho-II 
phase at room temperature scanned along h (left), k (mid- 
dle) and I (right) in crystal #1. The thermal treatment of 
the crystal corresponds to state 2 as described in the text. 
The lines are least square fits to equation (|l|). The peak at 
h = 2.68 is an Al (2 2 0) powder-line. 



FIG. 2. Half width half maximum (HWHM) of the ortho-II 
superlattice peak along h, k and I plotted as function of sam- 
ple mosaicity in the ac-plane of four differently grown crystals. 
The growth techniques are described in the text. 



FIG. 3. Peak to background ratio of the ortho-II superlat- 
tice reflection (2.5 5) plotted as function of the peak width 
measured at crystal #1 with x = 0.50 after the different ther- 
mal treatments (labeled 1 ... 6) that are described in the text. 
The lines are fits to a quadratic dependence of the intensity 
on the width. 



FIG. 4. Top: (2.5 5) peak intensity measured by tj-scans 
for crystal #1 with x — 0.50 after annealing at 80 °C for 70 
days and cooling to room tem perat ure by 1 °C/hour (anneal- 
ing procedure 1 in Subsection [II A). Open symbols mark the 



data obtained during heating, closed symbols those from cool- 
ing. The inset shows the slope of the data for heating (solid 
line) and those for cooling (dashed line)normalized to unity. 
Middle: Half width half maximum of the superstructure re- 
flection. The line is a fit to the critical behavior of a 3-di- 
mensional Ising model. Bottom: Exponent of the Lorentzian 
scattering function. The exponent of the data above 95°C are 
fixed to y' = 1. 



FIG. 5. Ortho-Ill superlattice reflections measured on a 
crystal with oxygen composition x — 0.72 at room tempera- 
ture. Left: (h 5)-scan, middle: (8/3 k 0)-scan and right: 
scan along (7/3 I) (open circles) and (8/3 I) (filled cir- 
cles). The lines shown with the h- and fc-scan are fits to a 
Lorentzian function (y 1 = 1). 



FIG. 6. Ortho-Ill transition at the oxygen composition 
x = 0.72 where the peak intensity and the width of the 
(8/3 5)-reflection were measured by /i-scans. The open 
circles show the heating data and the filled circles the data 
during cooling. The line in the lower graph shows the behav- 
ior described by the 2D Ising model. 



FIG. 7. Scans along h (left), k (middle) and I (right) at 
room temperature for a crystal prepared to x — 0.62. The 
peaks in the (h 0)-scan at h = 2.4 and h — 2.6 result from 
the ortho-V phase. The scattering signal at h = 2.5 indicates 
the presence of the ortho-II phase. The weak peak at h — 2.23 
is an Al-powder line originating from the sample holder and 
the bump at h = 2.83 is unidentified. At h = 2.1 and h — 2.9 
tails of the fundamental Bragg reflections are observed. In 
the middle part of the figure white circles show a scan along 
(2.4 k 0) and black squares a (2.5 k 0)-scan. The right hand 
side shows a (2.5 Z)-scan. The arrows indicate the position 
of Al-powder lines. 



FIG. 8. Temperature dependence of the mixed phase of 
ortho-II and ortho-V at the oxygen composition of x = 0.62. 
The top part plots the peak intensity of the ortho-II reflections 
(circles) and the ortho-V reflections (squares) . Open and gray 
shaded symbols mark the data obtained during heating, the 
black filled symbols the cooling data. The bottom part of the 
figure shows the half width at half maximum (HWHM) of the 
reflections of the two phases along h. 
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FIG. 9. Scan along (h 0) (left), (2.63 k 0) (middle) and 
(2.63 I) (right) in a crystal with oxygen concentration of 
x=0.67 at room temperature. The /i-scan on the left hand side 
shows peaks at approximately h — 2| and h = 2| indicating 
the ortho-VIII phase. The peak at 2.23 is an Al-powder line. 
Lines shown together with the h- and fc-scans are Lorentzian 
fits. The variation of the diffracted intensity along I is only 
weakly q-dependent and similar to that of the ortho-Ill phase. 



FIG. 10. Temperature dependence of the ortho-VIII phase 
at x = 0.67 measured by /i-scans on the (2| 3) superstruc- 
ture reflection. The top figure shows the peak intensity, the 
middle part the half width half maximum together with a fit 
to a critical behavior according to equation H with v = 0.79(3) 
and the bottom the peak position. The shift in the peak po- 
sition indicates a transformation from ortho-VIII to ortho-V 
structure between 50°C and 70°C. At 150°C the peak position 
corresponds to the ortho-Ill structure. 



FIG. 11. The structural phase diagram of YBCO. The 
structural phases and their transition temperatures are la- 
beled: teragonal (T), ortho-I (OI,D), ortho-II (OII,«), or- 
tho-Ill (OIII,o), ortho-V (OV,A) and ortho-VIII (OVIII,o). 
Solid lines are guides to the eye. The dashed lines are pre- 
dictions from the ASYNNNI model. Th* data for the T-OI 
transition (D) are from Andersen et al.a. The Ton transi- 
tion temperatures for x — 0.35 and x— 0.36 are from Poulsen 
eipaZ.tj, the upper data set for 1=0.50 is from Schleger et 
alu, and the TrSiiii transition temperature at a;=0.77 is from 
Schleger et al.c3. 



FIG. 12. Top: Half width half maximum for the indicated 
structures as function of oxygen content. Bottom: peak to 
background ratio. All lines are guides to the eye. 
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